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1 Introduction 



Intermediate energy heavy-ion collisions produce rich amount of information on the cor- 
relations and fluctuations and eventually on the dynamics and interactions among the 
nucleons. The breaking of nuclei, i.e., multifragmentation, is one of the rare phenomena 
that has attracted major attention in recent years [1]. The physics behind multifragmen- 
tation is so complicated that many different theoretical approaches have been developed 
[1-4]. Since no theoretical model simulates fragments, one needs afterburners to iden- 
tify clusters. Since correlations and fluctuations are the main features of the molecular 
dynamics model, the quantum molecular dynamics (QMD) model is very successful in 
explaining the phenomena of multifragmentation. Once the phase space is accessible, 
one generally clusterize the phase space with simple spatial correlation method where 
one binds the nucleons in a fragment they lie within a distance of 4 fm. This method is 
known as minimum spanning tree (MST) method [5]. At the same time fragments formed 
in MST method will be highly unstable (especially in central collisions) as there the two 
nucleons may not be well formed and therefore can be unstable that will decay after a 
while. In order to filter out such unstable fragments, we impose another cut in terms of 
relative momentum of nucleons. This method, dubbed as minimum spanning tree with 
momentum cut (MSTP) method was discussed by Puri et al. [6]. In our recent work, we 
study the role of momentum cut on fragment structure [7]. We also studied the role of 
colliding geometry on the fragmentation when momentum cut is being imposed. No study 
exits in literature to see the role of momentum cut on the various fragment properties 
like the rapidity distribution, p t spectra and E rat . So the the present paper, we plan 
to see the role of momentum cut on various fragment properties and to investigate how 
these properties vary with impact parameter. The present study is carried out within the 
framework of QMD model [2, 3] which is described in the following section. 
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2 The Formalism 



2.1 Quantum Molecular dynamics (QMD) model 

We describe the time evolution of a heavy-ion reaction within the framework of Quantum 
Molecular Dynamics (QMD) model [2, 3] which is based on a molecular dynamics picture. 
This model has been successful in explaining collective flow [8], elliptic flow [9], multifrag- 
mentation [10] as well as dense and hot matter [11]. Here each nucleon is represented by 
a coherent state of the form 

3 

<t>a{xi,t) = (J^J * e -(*i-*«(*))V p «( Xl -*»>e-^. (1) 

Thus, the wave function has two time dependent parameters x a and p a . The total n-body 
wave function is assumed to be a direct product of coherent states: 

= <f>a(xi,X a ,p a ,t)<j>/3(x 2 ,Xi3,pp,t)...., (2) 

where antisymmetrization is neglected. One should, however, keep in the mind that the 
Pauli principle, which is very important at low incident energies, has been taken into 
account. The initial values of the parameters are chosen in a way that the ensemble 
(A T +A P ) nucleons give a proper density distribution as well as a proper momentum 
distribution of the projectile and target nuclei. The time evolution of the system is 
calculated using the generalized variational principle. We start out from the action 

S= / <2 £[0,0*]rfr, (3) 



ti 



with the Lagrange functional 

C 



dt 



(4) 



where the total time derivative includes the derivatives with respect to the parameters. 
The time evolution is obtained by the requirement that the action is stationary under the 
allowed variation of the wave function 

5S = 5 [ * £[<f),<i)*}dt = 0. (5) 

If the true solution of the Schrodinger equation is contained in the restricted set of 
wave function <p a (xi,x a ,p a ) , this variation of the action gives the exact solution of the 
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Schrodinger equation. If the parameter space is too restricted, we obtain that wave func- 
tion in the restricted parameter space which comes close to the solution of the Schrodinger 
equation. Performing the variation with the test wave function (2), we obtain for each 
parameter A an Euler-Lagrange equation; 

d dC dC n (r\ 

dt^;-d\ =0 - (6) 

For each coherent state and a Hamiltonian of the form, 

H = J2 a r ^a J r\ K^J , the Lagrangian and the Euler-Lagrange function can be 
easily calculated 
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± a = ^ + V Pa J2(V a p), (8) 

m p 

p. = -Vi.E(u ( 9 ) 

Thus, the variational approach has reduced the n-body Schrodinger equation to a set of 6n- 
different equations for the parameters which can be solved numerically. If one inspects the 
formalism carefully, one finds that the interaction potential which is actually the Bruckner 
G-matrix can be divided into two parts: (i) a real part and (ii) an imaginary part. The 
real part of the potential acts like a potential whereas imaginary part is proportional to 
the cross section. 

In the present model, interaction potential comprises of the following terms: 

v af) = vL + vL + Vcoui + v Yuk (io) 

Vioc is the Skyrme force whereas Vcoui, Vy u k and Vmdi define, respectively, the Coulomb, 
and Yukawa potentials. The Yukawa term separates the surface which also plays the role 
in low energy processes like fusion and cluster radioactivity [12]. The expectation value 
of these potentials is calculated as 

VL = J / a (p a ,r a ,t)/ /3 (p /3 ,r /3 ,t)\/ / (2) (r a ,r /3 ) 

XCpTaCpTpcPpacPpp, (11) 
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Vile = J /a(p a ,r a ,t)/ /3 (p / 3,r / 3,t)/ 7 (p 7 ,r 7 ,t) 
xV} 3 \r a , Tp, r 7 )rf 3 r a rf 3 r /3 rf 3 r 7 

xd 3 p a d 3 P/3 d 3 p 7 . (12) 

where f a (p a , t) is the Wigner density which corresponds to the wave functions (eq. 2). 
If we deal with the local Skyrme force only we get 



A T +A P 
~Y Skyrme ^ ^ 

a=l 



(13) 



Here A, B and C are the Skyrme parameters which are defined according to the ground 
state properties of a nucleus. Different values of C lead to different equations of state. A 
larger value of C (= 380 MeV) is often dubbed as stiff equation of state. The finite range 
Yukawa (Vyuk) and effective Coulomb potential (Vcw) read as: 

exp{-\ri - ijl}//! 



v ™-2> | ri _ r .| /Ai > ( 14 ) 



^ ri - rj 

The Yukawa interaction (with t 3 = -6.66 MeV and pL = 1.5 fm) is essential for the 
surface effects. The relativistic effect does not play role in low incident energy of present 
interest [13]. 

The phase space of nucleons is stored at several time steps. The QMD model does 
not give any information about the fragments observed at the final stage of the reac- 
tion. In order to construct the fragments, one needs clusterization algorithms. We shall 
concentrate here on the MST and MSTP methods. 

According to MST method [5], two nucleons are allowed to share the same fragment 
if their centroids are closer than a distance r min , 

\ri - rj| < r min . (16) 

where ri and rj are the spatial positions of both nucleons and r min taken to be 4fm. 

For MSTP method,we impose a additional cut in the momentum space, i.e., we allow 
only those nucleons to form a fragment which in addition to equation(16) also satisfy 

|Pi-Pj| <Pmin, (17) 

where p min = 150 MeV/c. 
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Figure 1: The rapidity distribution for free nucleons and LCPs for the reaction of 12 C+ 12 C 
at incident energy of 100 MeV/nucleon at central (left panel) and peripheral (right) ge- 
ometries with MST and MSTP methods. 

3 Results and Discussion 

We simulated the reactions of 12 C+ 12 C and 40 Ca+ 40 Ca at 100 MeV/nucleon at central 
and peripheral colliding geometries, i.e., at b = 0.0 and 0.8, respectively. We use a soft 
equation of state with standard energy- dependent Cugon cross section. 
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Figure 2: Same as Fig. 1 but for the reaction of Ca+ Ca. 

In figure 1, we display the rapidity distribution of free nucleon and LCP's for the 
reaction of 12 C+ 12 C at energy of 100 MeV at central (left panel) and peripheral (right) 
colliding geometry. The solid and dashed lines indicate the calculations of MST and 
MSTP methods, respectively. From the figure, we see that there is quantitative difference 
in the results of MST and MSTP methods, though qualitatively, both methods give similar 
behaviour of the rapidity distribution of nucleons and fragments. 

For central collisions (left panel), we see that the peak of dN/dY plot is pronounced 
for the MSTP method, thus, indicating enhanced production of free nucleons in MSTP 
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Figure 3: The (l/(MeV/c) 2 ) as a function of transverse energy p t for the free nucle- 
ons and LCPs for the reaction of 12 C+ 12 C at central (left panel) and peripheral (right) 
geometries with MST and MSTP methods. Lines have same meaning as in Fig. 1. 

method as compared to MST method. This is due to the fact that in the MST method, 
we have a single big fragment because of no restriction is being imposed on the relative 
momentum of nucleons forming fragments. The production of LCP's is more with MSTP 
method compared to MST method which is supported by Ref. [6,7]. At peripheral 
collisions, the behaviour of the rapidity plots of free nucleons is similar as for the central 
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Figure 4: Same as Fig. 3 but for the reaction of Ca+ Ca. 

one whereas the trend reverses for the LCP's plot. For peripheral collisions, we have 
greater production with MSTP method. 

In figure 2, we display the rapidity distribution of free nucleons, LCP's and IMF's for 
the reaction of 40 Ca+ 40 Ca. Left(right) panels display the results for b/b maa; =0.0 (0.8). We 
find similar behaviour of free nucleons and LCP's as reported for the reaction of 12 C+ 12 C. 
The IMF's also follow the similar trend as for LCP's, i.e., we have more (less) production 
of IMF's with MST method at central (peripheral) collisions. 

In figures 3 and 4, we display dN/pjdpt versus p t for the reaction of 12 C+ 12 C and 
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Figure 5: The time evolution of E rat for free nucleons and LCPs for the reaction of 
12 C+ 12 C. 

40 Ca+ Ca, respectively. We see that dN/p t dpj spectra follow the similar behaviour for 
both MST and MSTP methods. We have a higher peak in the spectra of free nucleons 
with MST method at both the colliding geometries. The difference between MST and 
MSTP methods in spectra of LCP's is less significant. Similar behaviour is also observed 
for the reaction of 40 Ca+ 40 Ca. 

In figure 5, we display the time evolution of E rat of free nucleons and LCP's for 
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Figure 6: Same as Fig. 5 but for the reaction of 40 Ca+ 40 Ca. 

the reaction of 12 C+ 12 C at central (left panel) and peripheral (right) colliding geometry. 
For MST and MSTP methods, we find a significant difference between MST and MSTP 
methods for both free nucleons and LCP's. The difference is more for the central collisions 
as compared to the peripheral one. 

In figure 6, we display time evolution of E rat of free nucleons, LCP's and IMF's for 
the reaction of 40 Ca+ 40 Ca at central (left panel) and peripheral (right) collisions. The 
solid (dashed) lines represent the results of MST (MSTP) method. From figure, we find 
a significant difference of E rat with MST and MSTP method as in case of the reaction of 
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12 C+ 12 C. We also find that the difference between MST and MSTP reduces at peripheral 
colliding geometries. 



4 Summary 

Using the quantum molecular dynamic model, we studied the role of momentum corre- 
lations in the properties of fragments. This was achieved by imposing cut in momentum 
space during the process of clusterization. We find that this cut yields significant differ- 
ence in the fragment properties of system at all colliding geometries. 
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